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Abstract 

, In this article, we perform a systematic study of the mass spectrum of the scalar 

(3JT)' doubly charmed and doubly bottom tetraquark states using the QCD sum rules. 
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r—i 1 Introduction 

The Z(4430) and the Z(4050), Z(4250) observed in the decay modes ^'tt + and xa^ + 
respectively by the Belle collaboration are the most interesting subjects [H [2j [3]. We 
can distinguish the multiquark states from the hybrids or charmonia with the criterion of 
non-zero charge. They can't be pure cc states due to the positive charge, and must be 

Cf"} ■ some special combinations of the ccud tetraquark states, irrespective of the molecule type 

and the diquark-antidiquark type. If those states are confirmed in the future, they are 

qq | excellent candidates for the heavy tetraquark states of the QqQq' type. It is interesting 

to explore the possibility that whether or not there exist doubly heavy tetraquark states 
of the QQqq' type. 

On the other hand, the QCD sum rules is a powerful theoretical tool in studying 
the ground state hadrons [4j [5]. In the QCD sum rules, the operator product expansion 
t— I . is used to expand the time-ordered currents into a series of quark and gluon condensates 

which parameterize the long distance properties of the QCD vacuum. Based on the quark- 
hadron duality, we can obtain copious information about the hadronic parameters at the 
phenomenological side [H [5] . 

There have been several successful applications of the QCD sum rules in studying the 
hidden charmed and hidden bottom tetraquark states (QqQq' type). In Refs.jH [8], we 
study the mass spectrum of the scalar hidden charmed and hidden bottom tetraquark 
states in a systematic way using the QCD sum rules, and identify the Z(4250) tentatively 
as a scalar tetraquark state of the diquark-antidiquark type; while in Ref.[9] the Z(4050) 
and Z(4250) are interpreted as the D\D molecular state. In Refs.[T0~l ITT] , we study 
the mass spectrum of the vector hidden charmed and hidden bottom tetraquark states 
systematically using the QCD sum rules. In Ref.|12j. we perform a systematic study of 
the mass spectrum of the axial-vector hidden charmed and hidden bottom tetraquark 
states using the QCD sum rules, and identify the Z(4430) tentatively as an axial-vector 
tetraquark state of the diquark-antidiquark type. In Refs.[13] [H], Lee et al study the 
Z(4430) with the QCD sum rules and observe that the Z(4430) maybe a - molecular 
type or diquark-antidiquark type tetraquark state. In Ref.[15j, Chen et al study the 
hidden charmed and hidden bottom tetraquark states in details with the QCD sum rules. 
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In Ref.|16j. Navarra et al use the QCD sum rules to study the possible existence 
of the doubly heavy tetraquark states QQud with J p = 1 + . There have been several 
other theoretical approaches in studying the doubly heavy tetraquark states, such as the 
potential models and QCD inspired potential models [EJ HH1 [H] , solving the four-body 
problem within a non-relativistic quark model |20[ [21] , the variational method combined 
with a non-relativistic potential model [22] , the chiral constituent quark model |23[ [24] , the 
semi-empirical mass relations [25], the relativistic quark model based on a quasipotential 
approach in QCD |26j . etc. Whether or not there exist the doubly charmed or doubly 
bottom tetraquark configurations is of great importance itself, because it provides a new 
opportunity for a deeper understanding of the low energy QCD. 

It is interesting to study the mass spectrum of the doubly heavy tetraquark states 
(QQqq' type) with the QCD sum rules, and make an independent estimation from QCD. 
In Refs.pTl l28j [29] . we study the mass spectrum of the \ and | doubly heavy baryon 
states (QQq type) in a systematic way using the QCD sum rules. In this article, we extend 
our previous works to study the mass spectrum of the scalar doubly charmed and doubly 
bottom tetraquark states in a systematic way with the QCD sum rules. 

We take the diquarks as the basic constituents following Jaffe and Wilczek [30 S [3T], and 
construct the doubly heavy tetraquark states with the diquark and antidiquark pairs. The 
diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vector 07^75, axial- 
vector C7 M and tensor Ca^, where C is the charge conjunction matrix. The structures 
C7 M and Ca^ u are symmetric while the structures C75, C and 6*7^75 are antisymmetric. 

The scattering amplitude for one-gluon exchange in an SU{N C ) gauge theory is pro- 
portional to 

T ki T ij = - -^-^ (SjhSil ~ S ik Sji) + ° m 1 (S jk 6u + SikSji) , (1) 

where the T a is the generator of the gauge group, and the i, j and k, I are the color indexes 
of the two quarks in the incoming and outgoing channels respectively. For N c = 3, the 
negative sign in front of the antisymmetric antitriplet indicates the interaction is attractive, 
while the positive sign in front of the symmetric sextet indicates the interaction is repulsive 
[32j. On the other hand, the scattering amplitude for one-gluon exchange in the Dirac 
spinor space is proportional to 

{l^ijWkl = -h5C)ik{Cj 5 )lj + {C) ik {C)ij + i(757aC)ifc(C7 Q 75)ii 

— (laCMC^u , (2) 

the negative sign in front of the scalar and axial-vector channels indicates the interaction 
is attractive. 

For the doubly heavy quark system with the same flavor Q t CTQ :) , where the T de- 
note the Dirac matrixes 1, 75, 7^, 7^75 and a^, the color indexes i and j should be 
antisymmetric, i.e. 

Q l CTQ j ~ e ijk Q i CTQ j . (3) 

In the case of the antisymmetric structures C75, C and 6*7^75, the fermi statistics forbids 
the formulation of the diquark states. 
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In this article, we use the symmetric structure Cj^ to construct the interpolating 
currents J(x) to study the doubly charmed and doubly bottom tetraquark states Z: 



Jqq{%) 
•Jqs (•£■) 
Jss 



e ij h imn Qj(x)C^Q k (x)q m (x)j»CqT(x), 
e^ k e mn Qj(x)C^Q k (x)q m (x)^Csl(x), 
e^ k e imn Q]{x)C ltl Q k {x)s m {x)^Cs T n {x) : 



(4) 



where q = u,d. In the isospin limit, the interpolating currents J{x) result in three distinct 
expressions for the correlation functions II (p), which are characterized by the number of 
the s quark they contain. In Refs.[71 [11], we observe that the ground state masses of the 
scalar and vector hidden charmed and hidden bottom tetraquarks are characterized by 
the number of the s quarks they contain, Mq < M s < M ss ; the energy gap between M$ 
and M ss is about (0.05 — 0.15) GeV. In this article, we study the interpolating currents 
which contains zero and two s quarks for simplicity. 

Lattice QCD calculations for the light flavors indicate that the strong attraction in the 
scalar diquark channels favors the formation of good diquarks, the weaker attraction (the 
quark-quark correlation is rather weak) in the axial-vector diquark channels maybe form 
bad diquarks, the energy gap between the axial- vector and scalar diquarks is about | of the 
A-nucleon mass splitting, i.e. ~ 0.2 GeV [331 134] , which is expected from the hypersplit- 
ting color-spin interaction m l m . Ti ■ Tj&i ■ aj [31J. The coupled rainbow Dyson-Schwinger 
equation and ladder Bethe-Salpeter equation also indicate such an energy hierarchy [35] . 
Comparing with the spin independent term Ti-Tj, the contribution from the hypersplitting 
color-spin interaction mm . Ti ■ Tjffi ■ aj is greatly suppressed by the inverse constituent 
quark masses. It is possible to form axial-vector diquark states, although the hypersplit- 
ting color-spin interaction m 1 m Ti-Tj&i-aj is repulsive in this channel. If we take the scalar 
light diquark states as the basic constituents, additional relative P- waves are needed to 
obtain the correct zero spin. In the conventional quark models, additional P-wave excita- 
tion costs about 0.5 GeV, the ground states should be constructed with the axial-vector 
antidiquark states. 

The article is arranged as follows: we derive the QCD sum rules for the scalar doubly 
charmed and doubly bottom tetraquark states Z in Sect. 2; in Sect. 3, we present the 
numerical results and discussions; and Sect. 4 is reserved for our conclusions. 

2 QCD sum rules for the scalar tetraquark states Z 

In the following, we write down the two-point correlation functions H(p) in the QCD sum 
rules, 



where the J(x) denotes the interpolating currents J gg (x) and J ss (x). 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J(x) into the correlation functions II(p) to obtain the 
hadronic representation [HE]. After isolating the ground state contribution from the pole 
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term of the Z, we get the following result, 
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n(p) = jit^ + -- W 

where the pole residue (or coupling) Xz is defined by 

X z = (0|J(0)|Z(p)>. (7) 

After performing the standard procedure of the QCD sum rules, we obtain the following 
two sum rules: 

A|e m 2 " = / dsp z {s)e TP , (8) 

the explicit expressions of the spectral densities pz(s) are presented in the appendix, the 

is the continuum threshold parameter and the M 2 is the Borel parameter. We can 
obtain two sum rules in the ccqq and bbqq channels with a simple replacement m s — > m q , 
(ss) -> (qq) and (sg s aGs) -> (qg s oGq). 

We carry out the operator product expansion to the vacuum condensates adding up to 
dimension-10. In calculation, we take vacuum saturation for the high dimension vacuum 
condensates, they are always factorized to lower condensates with vacuum saturation in 
the QCD sum rules, factorization works well in large iV c limit. In reality, N c = 3, some 
ambiguities may come from the vacuum saturation assumption. 

We take into account the contributions from the quark condensates, mixed condensates, 
and neglect the contributions from the gluon condensate. The gluon condensate ( " s ) 
is of higher order in a s , and its contributions are suppressed by very large denominators 
and would not play any significant role for the light tetraquark states |36|, [37] , the heavy 
tetraquark state [6] and the heavy molecular states [38| I3U]. 

In the special case of the y(4660) (as a ip'f b(980) bound state) and its pseudoscalar 
partner rj' c fo(980), the contributions from the gluon condensate ( as ) are rather large 
[4171 SI]- If we take a simple replacement s(x)s(x) — > (ss) and [u(x)u(x) + d(x)d(x)~j — > 
2(qq) in the interpolating currents, the standard heavy quark currents Q(x)^ lJi Q(x) and 
Q(x)i'j 5 Q(x) are obtained, where the gluon condensate ( as r ) plays an important rule 
in the QCD sum rules [4]. The interpolating currents constructed from the diquark- 
antidiquark pairs do not have such feature. 

We also neglect the terms proportional to the m u and m^, their contributions are of 
minor importance due to the small values of the u and d quark masses. 

Differentiating the Eq.(8) with respect to then eliminate the pole residues Xz, we 
can obtain the sum rules for the masses of the Z, 

M Z = JAZ f^ ^l . (9) 

j A z z dsp z (s)e ^ 

3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = —(0.24 ± 0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2)(qq), (qg s aGq) = m 2 (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 0.2) GeV 2 , 
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m s = (0.14 ± 0.01) GeV, m c = (1.35 ± 0.10) GeV and m b = (4.8 ± 0.1) GeV at the energy 
scale fi = lGeV [11 El S2]. 

The Q-quark masses appearing in the perturbative terms are usually taken to be 
the pole masses in the QCD sum rules, while the choice of the tbq in the leading-order 
coefficients of the higher-dimensional terms is arbitrary [A3\ |44J. The MS mass m c (m 2 ) 

relates with the pole mass m c through the relation m c (m 2 ) = rh c 1 + CFaa J- m ^ -{-■■■ 

In this article, we take the approximation m c (m 2 ) ~ m c without the a s corrections for 
consistency. The value listed in the Particle Data Group is m c (m 2 ) = 1.27^o'ii GeV [43] . 
it is reasonable to take rh c = m c (1 GeV 2 ) = (1.35 ± 0.10) GeV. For the b quark, the MS 
mass mb(m 2 ) = 4.20^q'q7 GeV [45], the gap between the energy scale fx = 4.2 GeV and 
1 GeV is rather large, the approximation mj, ~ mb(m 2 ) ~ mj,(l GeV 2 ) seems rather crude. 
It would be better to understand the quark masses m c and nib we take at the energy scale 
[i 2 = 1 GeV 2 as the effective quark masses (or just the mass parameters). Our previous 
works on the mass spectrum of the heavy and doubly heavy baryon states indicate such 
parameters can lead to satisfactory results [271 EH EHJ HH H7J . 

In calculation, we also neglect the contributions from the perturbative corrections. 
Those perturbative corrections can be taken into account in the leading logarithmic ap- 
proximations through anomalous dimension factors. After the Borel transform, the effects 
of those corrections are to multiply each term on the operator product expansion side by 

the factor, 



where the Tj is the anomalous dimension of the interpolat- 
ing current J(x) and the To n is the anomalous dimension of the local operator O n (0). 

We carry out the operator product expansion at a special energy scale \i 2 = 1 GeV 2 , and 

2r,-r„„ 



set the factor 



a 3 (M 2 ) 
~~aAp I ) 



~ 1, such an approximation maybe result in some scale 
dependence and weaken the prediction ability. In this article, we study the scalar doubly 
charmed and doubly bottom tetraquark states systemically, the predictions are still robust 
as we take the analogous criteria in those sum rules. 

In the conventional QCD sum rules [HE], there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter M 2 and 
threshold parameter sq. We impose the two criteria on the scalar doubly heavy tetraquark 
states to choose the Borel parameter M 2 and threshold parameter s$. 

The vacuum condensates of the high dimension play an important role in choosing the 
Borel parameter M 2 . The condensate of the highest dimension (sg s aGs) 2 is counted as 

O(^p-), C(ir#) or O(j^), and the corresponding contributions are greatly enhanced at 
small M 2 , and result in rather bad convergent behavior in the operator product expansion, 
we have to choose large Borel parameter M 2 . We insist on taking into account the high 
dimensional vacuum condensates, as the interpolating current consists of a (heavy )diquar k- 
( light )antidiquark pair, one of the highest dimensional vacuum condensates is {ss) 2 ( as ^ G ), 
we have to take into account the condensate (sg s aGs) 2 for consistence. 

The contributions from the high dimension vacuum condensates in the operator prod- 
uct expansion are shown in Fig.l, where (and thereafter) we use the (ss) to denote the 
quark condensates (qq), (ss) and the (sg s aGs) to denote the mixed condensates (qg s aGq), 
(sg s aGs). From the figures, we can see that the contributions from the high dimension 
condensates are very large and change quickly with variation of the Borel parameter at 
the values M 2 < 2.6 GeV 2 and M 2 < 7.2 GeV 2 in the doubly charmed and doubly bot- 
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torn channels respectively, such an unstable behavior cannot lead to stable sum rules, our 
numerical results confirm this conjecture, see Fig. 3. 

At the values M 2 > 2.6 GeV 2 and s > 25GeV 2 , 24GeV 2 , the contributions from the 
(ss) 2 + (ss)(sg s aGs) term are less than 9%, 23% in the channels cess, ccqq respectively; 
the contributions from the vacuum condensate of the highest dimension (sg s oGs) 2 are 
less than 4%, 5% in the channels cess, ccqq respectively; we expect the operator product 
expansion is convergent in the doubly charmed channels. 

At the values M 2 > 7.2 GeV 2 and s > 140 GeV 2 , 138 GeV 2 , the contributions from the 
(ss) 2 + (ss)(sg s aGs) term are less than 6%, 16% in the channels bbss, bbqq respectively; 
the contributions from the vacuum condensate of the highest dimension (sg s oGs) 2 are 
less than 5%, 8% in the channels bbss, bbqq respectively; we expect the operator product 
expansion is convergent in the doubly bottom channels. 

In this article, we take the uniform Borel parameter M^ in , i.e. M^ in > 2.6 GeV 2 and 
^min — 7-2 GeV 2 in the doubly charmed and doubly bottom channels respectively. 

In Fig. 2, we show the contributions from the pole terms with variation of the Borel 
parameters and the threshold parameters. The pole contributions are larger than (or 
equal) 50%, 47% at the value M 2 < 3.3 GeV 2 and s > 25 GeV 2 , 24 GeV 2 in the channels 
cess, ccqq respectively, and larger than (or equal) 52%, 50% at the value M 2 < 8.2 GeV 2 
and so > 140 GeV 2 , 138 GeV 2 in the channels bbss, bbqq respectively. Again we take the 
uniform Borel parameter M^ ax , i.e. M^ ax < 3.3 GeV 2 and M^ ax < 8.2 GeV 2 in the 
doubly charmed and doubly bottom channels respectively. 

In this article, the threshold parameters are taken as so = (26±l)GeV 2 , (25±l)GeV 2 , 
(142 ± 2) GeV 2 , (140 ± 2) GeV 2 in the channels cess, ccqq, bbss, bbqq respectively; the 
Borel parameters are taken as M 2 = (2.6 - 3.3) GeV 2 and (7.2 - 8.2) GeV 2 in the doubly 
charmed and doubly bottom channels respectively. In those regions, the pole contributions 
are about (50 - 80)%, (47 - 78)%, (52 - 71)%, (50 - 70)% in the channels cess, ccqq, bbss, 
bbqq respectively; the two criteria of the QCD sum rules are fully satisfied [4, 5j. 

From Fig.2, we can see that the Borel windows M^ ax — M^ in change with variations of 
the threshold parameters so- In this article, the Borel windows are taken as 0.7 GeV 2 and 
1.0 GeV 2 in the doubly charmed and doubly bottom channels respectively; they are small 
enough. If we take larger threshold parameters, the Borel windows are larger and the 
resulting masses are larger, see Fig. 3. In this article, we intend to calculate the possibly 
lowest masses which are supposed to be the ground state masses by imposing the two 
criteria of the QCD sum rules. 

Taking into account all uncertainties of the relevant parameters, finally we obtain the 
values of the masses and pole resides of the scalar doubly heavy tetraquark states Z, which 
are shown in Figs. 4-5 and Table 1. 

In this article, we calculate the uncertainties 5 with the formula 




where the / denote the hadron mass Mz and the pole residue Xz , the x% denote the relevant 
parameters m c , m&, (qq), (ss), ■ ■ ■ . As the partial derivatives J^- are difficult to carry out 

analytically, we take the approximation (j^J (xi — Xi) 2 ~ [f(xi ± Axi) — f(xi)] 2 in the 
numerical calculations. 
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Figure 1: The contributions from different terms with variation of the Borel parameter 
M 2 in the operator product expansion. The (I) and (II) denote the contributions from 
the (ss) 2 + (ss}(sg s aGs) term and the {sg s aGs} 2 term respectively. The A, B, C and D 
denote the channels cess, ccqq, bbss and bbqq respectively. The notations a, /3, 7, A, p 
and r correspond to the threshold parameters s = 22 GeV 2 , 23 GeV 2 , 24 GeV 2 , 25 GeV 2 , 
26 GeV 2 and 27 GeV 2 respectively in the doubly charmed channels; while they correspond 
to the threshold parameters s = 134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 140 GeV 2 , 142 GeV 2 
and 144 GeV 2 respectively in the doubly bottom channels. 
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Figure 2: The contributions of the pole terms with variation of the Borel parameter 
M 2 . The A, B, C and D denote the channels cess, ccqq, bbss and bbqq respectively. 
The notations a, (3, 7, A, p and r correspond to the threshold parameters sq = 22 GeV 2 , 
23 GeV 2 , 24 GeV 2 , 25 GeV 2 , 26 GeV 2 and 27 GeV 2 respectively in the doubly charmed 
channels; while they correspond to the threshold parameters so = 134 GeV 2 , 136 GeV 2 , 
138 GeV 2 , 140 GeV 2 , 142 GeV 2 and 144 GeV 2 respectively in the doubly bottom channels. 
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bbss 
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0.825 ±0.180 


bbqq 
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0.660 ±0.164 



Table 1: The masses and the pole residues of the scalar doubly heavy tetraquark states. 
The masses are in unit of GeV and the pole residues are in unit of GeV 5 . 
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Figure 3: The masses of the scalar doubly heavy tetraquark states with variation of 
the Borel parameter M 2 . The A, B, C and D denote the channels cess, ccqq, bbss 
and bbqq respectively. The notations a, f3, 7, A, p and r correspond to the threshold 
parameters s = 22 GeV 2 , 23 GeV 2 , 24 GeV 2 , 25 GeV 2 , 26 GeV 2 and 27 GeV 2 respectively 
in the doubly charmed channels; while they correspond to the threshold parameters sq = 
134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 140 GeV 2 , 142 GeV 2 and 144 GeV 2 respectively in the 
doubly bottom channels. 
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Table 2: The masses of the scalar doubly heavy tetraquark states, the star denotes the 
corresponding QqQq' type tetraquark states. The masses are in unit of GeV. 
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Figure 4: The masses of the scalar doubly heavy tetraquark states with variation of the 
Borel parameter M 2 . The A, B, C and D denote the channels cess, ccqq, bbss and bbqq 
respectively. 
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Figure 5: The pole residues of the scalar doubly heavy tetraquark states with variation 
of the Borel parameter M 2 . The A, B, C and D denote the channels cess, ccqq, bbss and 
bbqq respectively. 



11 



Naively, we expect the QQqq' and QqQq' type tetraquark states have degenerate masses 
as the color interactions are flavor blinded. However, we cannot obtain a relation between 
the corresponding interpolating currents by Fierz reordering in the color and Dirac spinor 
spaces as different flavor structures are concerned. Furthermore, additional contributions 
from the instanton configurations make the situation more complicated |48j . In Table 
2, we also present the values of the QQqq' and QqQq' type tetraquark states from the 
relativistic quark model based on a quasipotential approach in QCD [261 SHI EQ], the 
constituent diquark model plus the spin-spin interactions [HJ ES [53] , and the QCD sum 
rules [8]. The QCD sum rules indicate the QQqq' and QqQq' type tetraquark states have 
almost degenerate masses, while the central values of present predictions are larger than 
the corresponding ones from other theoretical models about (0.2 — 0.7) GeV. 

In Refs.[26j 091 [50], Ebert et al take the diquarks as bound states of the two quarks in 
the color antitriplet channel, and calculate their mass spectrum using a Schrodinger type 
equation, then take the masses of the diquarks as the basic input parameters, and study the 
mass spectrum of the heavy tetraquark states as bound states of the diquark-antidiquark 
system. In Refs.[5TJ [52], [53], Maiani et al take the diquarks as the basic constituents, 
examine the rich spectrum of the diquark-antidiquark states with the constituent diquark 
masses and the spin-spin interactions, and try to accommodate some of the newly observed 
charmonium-like and bottonium-like resonances not fitting pure cc and bb assignment. 
The predictions depend heavily on the assumption that the light scalar mesons ao(980) 
and /o(980) are tetraquark states, the basic parameters (constituent diquark masses) are 
estimated thereafter. In the conventional quark models, the constituent quark masses 
are taken as the basic input parameters, and fitted to reproduce the mass spectra of the 
well known mesons and baryons. However, the present experimental knowledge about 
the phenomenological hadronic spectral densities of the tetraquark states is rather vague, 
the constituent diquark masses and thereafter predictions cannot be confronted with the 
experimental data. 

If kinematically allowed, the scalar doubly heavy tetraquark states Z can decay to the 
heavy meson pairs with the Okubo-Zweig-Iizuka super-allowed "fall-apart" mechanism, 
i.e. QQqq ->■ QqQq and QQss -> QsQs. The thresholds for the DD, D S D S , D*D*, 
D*D*, BB, B S B S , B*B* and B*B* are about 3.74 GeV, 3.94 GeV, 4.02 GeV, 4.22 GeV, 
10.56 GeV, 10.73 GeV, 10.65 GeV and 10.83 GeV, respectively 05]. From Table 1, we can 
see that the strong decays Z cc qq — > DD, D*D* and Z ccSS — > D S D S , D*D* are kinematically 
allowed, the phase spaces are rather large; while the corresponding decays for the doubly 
bottom tetraquark states are kinematically forbidden. 

The doubly heavy tetraquark states can also decay to the baryon pairs with creation 
of the qq or ss pairs from the QCD vacuum, QQqq — > QQq' + q'qq, QQs + sqq, QQss — > 
QQq + qss, QQs + sss. However, the strong decays to 'EqqP, HqqA, VLqqY>, Q,qqY<*, 
EqqE', EqqE*, QqqCI and QqqQ* are kinematically forbidden or greatly suppressed. 
The scalar doubly charmed tetraquark states maybe have large widths, while the scalar 
doubly bottom tetraquark states maybe have very small widths. 

In 2002, the SELEX collaboration reported the first observation of a signal for the 
doubly charm baryon state H+ in the charged decay mode — >• A+K~ir + [51], and 
confirmed later by the same collaboration in the decay mode — > pD + K~ with measured 
mass Ms = 3518.9 ±0.9 Me V |55| . No other doubly heavy baryon states are observed. We 
use the masses of the \ and | doubly heavy baryon states Hqq, Sqq, Oqq, Qqq from 
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the QCD sum rules [271 [281 [29] . 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. In 
proton-proton collisions at y/s = 14TeV, the bb cross section is expected to be ~ 500/^6 
producing 10 12 bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 x 10 32 cm _2 sec _1 |56j . The scalar doubly heavy tetraquark states predicted in the 
present work may be observed at the LHCb, if they exist indeed. 

4 Conclusion 

In this article, we study the mass spectrum of the scalar doubly charmed and doubly bot- 
tom tetraquark states with the QCD sum rules in a systematic way. The mass spectrum 
are calculated by imposing the two criteria (pole dominance and convergence of the opera- 
tor product expansion) of the QCD sum rules. The present predictions can be confronted 
with the experimental data in the future at the LHCb. 

Appendix 

The spectral densities pz(s) at the level of the quark-gluon degrees of freedom: 
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